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If X is a stable process of index a G (0, 2) whose Levy measure 
has density cx~"~^ on (0, oo), and Si = supg^^^j^ Xt, it is known that 
P{Si >x)^ Aa-^x-" as a: ^ oo and P{Si <x) - Ba'^p-^x"" as 
X irO. [Here p — P{Xi > 0) and A and B are known constants.] It is 
also known that Si has a continuous density, m say. The main point 
of this note is to show that m{x) ^ Aa;"'-'^^^' as a- — oo and m{x) ^ 
Bx"''~^ as s^O. Similar results are obtained for related densities. 

1. Introduction and results. This paper was motivated by the following 
question which arises in connection with some problems in optimal stopping. 
Let X be a strictly stable process of index a G (0, 2) which has positive jumps 
so that its Levy measure has density 

, X /c+3;-(°+^), x>0, 

'^(") = lc!|xr("+^ x<o, 

where c+ > 0,c_ > 0. Assume also that X is not a subordinator. Let r^, for 
X > 0, denote the first passage time of X above level x, namely, 

Tx = ini{t:Xt > x}. 

Then does Tx have a density function? If so, how does it behave at zero and 
infinity? 

We prefer to rephrase this in terms of the maximum process defined by 
St = supg^t Xs which we can do because of the obvious identity 

P{tx >t) = P{St <x) = P{Si < xt-1), 

where rj = 1/a, and we have used the scaling property. It is easy to check 
that St has a continuous density function, mt say, and by scaling we have 
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mt{x) = t '^m{x/t'^) where m stands for rn-i, so we conclude that Tx has a 
continuous density function hx{t) given by 

(2) hx{t) = r]xt-i-^m{xt-'^). 

Moreover, we can read off the asymptotic behavior of at infinity and zero 
from the asymptotic behavior of m at zero and infinity, respectively. It is not 
difficult to make a conjecture about what this behavior is since it is known 
from [7], Theorem 4A (see also [5], Proposition 4, page 221) that 

(3) >x) >x) - asx^oo, 

a 

and from [7], Theorem 3A (see also Proposition 2, page 219 in [5]) that 

(4) P(Si < x) - —x'^P as X ; 0. 

ap 

Here p = P{Xi > 0), and A and B are explicitly known constants which can 
be expressed in terms of c+ and c„, and the result (4) is also valid in 
the spectrally negative case when c+ = < c_, a € (1, 2) and ap = l. [Recall 
that (3) does not hold in this case since P{Xi > x) is exponentially small at 
oo.] So the obvious conjecture is that, in these cases, 

(5) m(x) ^ ^x"*-""^"^-* as X — > oo 
and 

(6) m{x) Bx°f-^ asxiO. 

This question turns out to be closely related to a similar question about the 
density function p{x), say of 6^^^ where {9i^\o < s <t) denotes the stable 
meander of length t; informally, this is the stable process conditioned to stay 
positive up to time t. The following, which is the main result in this paper, 
confirms that these conjectures are true. 

Theorem 1. For any strictly stable process X which is such that \X\ is 
not a subordinator, (6) holds and there is a constant C € (0, cx)) such that 

(7) p{x) - Cx'^P as X i 0. 
If X also has positive jumps, (5) holds and also 

(8) p(x) ^ j4/3~^x"("+^) as x^^oo. 

Remark 2. Let /(x) = /i(x) where /j(x) denotes the density function 
of Xt] this exists, and is continuous and bounded on M because the char- 
acteristic function of Xt is in Li. A great deal is known about / (see, e.g., 
[13], pages 87-90). In particular, it is known that / has exactly the same 
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behavior as x — t- oo as that stated for m in (5) which is similar also to that 
of p in (8). However, its behavior as x J, is quite different to those of m 
and p, and, in fact, 

(9) lim/(2;) = L»e (0,oo), 

where the explicit value of D can be read off from (14.30) and (14.33) in 
[13]. 

Remark 3. It should be noted that some of our results are already 
known in the special case that a € (1,2) and c_ = 0, that is, the spectrally 
positive case. Here in [6] the semi-explicit form of the Wiener-Hopf factor- 
ization was exploited to show that m(x) = Yl'^ o„x°"~^ where the a.„ are 
given explicitly; since p = 1 — 1/a, this confirms (6) in this case. Another 
result in [6] is an expression for the Fourier transform of m, and this is used 
in [11] to show that in this case (5) also holds. 

Remark 4. If we knew in advance that m was ultimately monotone, as 
x — > cx) or as X — )■ 0, (5) and (6) would follow immediately from (3) and (4), 
but as we do not have this information, we have to use a different method. 
Also, in contrast to the case of X, we have no explicit knowledge of the 
characteristic function of S, so our method has to be rather indirect. 

Remark 5. We can read off from these results that the asymptotic 
behavior of the density of is given, for each fixed x G (0,oo), by 

h^{t) - riEx'^Pt-^P^^'^ as t ^ oo and 

(10) 

hx{t) ^ Arj/x"" asUO. 

2. Preliminaries on Levy processes. In our proofs we will use several 
identities for stable processes, all of which are in fact special cases of results 
valid for general Levy processes. It seems worthwhile to state these results 
in the general case. Therefore in this section alone, X will denote a generic 
Levy process which is not compound Poisson. S and / will be the associated 
supremum and infimum processes, we will write L and L* for the local times 
at zero of the reflected processes S — X and X — I, respectively, and n and 
n* will denote the characteristic measures of the excursions away from of 
these processes. We write e for a typical excursion, for its lifetime and 
^(t) = n*(C >t). Then 

P.(d.):=„-(e.gd.|C>t)= "''^'^t;'^^" 

7r*(t) 

is a probability distribution which, in the stable case, coincides with that of 
9[^\ the stable meander of length t, at time t. (See [5], page 234.) Our first 
result connects these quantities to the distribution of St- 
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Lemma 6. There is a constant, < k < oo, which depends only on the 
normalization of L and L* , such that, for t,x > 0, 



kP(S, > i) = [' 

(11) ^° 



n*{et-s > x,C > t — s){b + vr(s)) ds 



+ An*{et> x,C>t), 

where A denotes the drift, b the killing rate and vf(s) = n{(^ > s) the tail of 
the Levy measure of the increasing ladder time process T. 

Proof. First we note that P{St > x) = P{Xt -It> x) by duality. Next, 
if e^^'^ denotes the excursion of X — I starting at time r € G := {r : {X — I)r = 
0}, we see that 

P{Xt-It>x) = E{l^^,,_,,^^^^), 
where gt := sup(r <t:r ^ G}. But the compensation formula gives 

n*{et-s> x,Q>t- s)E{dL*{s)) 
and we conclude by showing that 

(12) kE{dL* {s)) = b ds + Tf {s) ds + A5o{ds). 

To see this, note that 



oo \ / fOO 



e{ I e-'i'dL*(s)]=Ei I e-«^*"Wdt' ^ 



where is the Laplace exponent of the decreasing ladder time process 
T* = L*~^. But, since the Levy measure of T is given by iT{ds) = n{C, S ds), 

q / e-^'inis) ds + bds + A6oids)) = (1 - e-^")n(C eds) + b + Aq 
Jo Jo 

the Laplace exponent of T. Then (12) follows by Laplace inversion and the 
identity ^*{q)^{q) = kq (see, e.g., [5], page 166). □ 

We also need two renewal- type equations which involve n*, but first we 
recall that if we write 

/•oo 

G{dt, dx) = / P{Ts € dt, Hs G dx) ds 
Jo 
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for the renewal measure in the increasing bivariate ladder process (T, H) of 
X, we have, for t > 0, a; > 0, 

(13) G{dt,dx) = k'n*{etedx,C>t)dt = k'^{t)Pt{dx)dt, 

where k' is a constant which depends only on the choice of the normalization 
of the local time L* of X — I. This fact comes from Theorem 5 of Alili and 
Chaumont [2], and the following result is also due to Alili and Chaumont. 

Lemma 7. For t,x>0 we have 

(14) tG{dt,dx)= I I G{du,dz)P^{Xt-uedx) 

Ju=0 Jz=0 

and 

(15) xG{dt,dx) = [ [ G{du,dz)^^^^P^{Xt-uedx). 

Ju=oJz=0 t — U 

The first of these statements is the obvious analogue of a result for random 
walks in [4], and, in fact, (14) is stated and proved in [1], the preliminary 
version of [2]. However it does not appear in [2] which does, however, contain 

(15) as (4.9). Unfortunately no proof is given there, but a proof can be found 
in [3]. 

3. Proof of the small time results. Throughout this section, X will be 
a stable process, so that the ladder time processes T and T* are stable 
subordinators of index p and (1 — p), respectively, and we will choose a nor- 
malization of the local times so that their exponents are given by ^{q) = 
and $*((/) = q^~^; this choice entails that n{( > t) = 7f(t) = i~^/r(l — p) and 
n*(C > t) = TT*{t) = tP-^/T{p). Note also that A = 6 = 0, and the constant k 
in Lemma 6 is equal to 1. Recall that rj = 1/a, and that the density ft of 
Xt is bounded and continuous on M, and has the scaling property 

(16) ft{x) = t-yixt-'^) where / = /i. 



Lemma 8. 9\ has a continuous and bounded density, pt{x) say, G{dt,dx) 
has a continuous bivariate density g{t,x), and these are related by 

(17) g{t,x) = k"tP-^pt{x) fort>0,x>0. 
Furthermore, m satisfies 

sinpvr ps{x) sinpvr s-^^pixs-i) 

(18) mix) = / — i — 7- :—ds = / — I — :—ds, 

^ ^ ^ ^ Jo s^-P{l-s)P vr Jo s^-Pil-s)P 

where we have written p for pi . 
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Proof. Note first that if we write ctq = inf {t : < 0} and introduce the 
measure, Qy{Xt S dx) := Py{Xt € dx, lt>0)= Py{Xt G dx, aQ> t), then, for 
x,y> 0, 



Py{Xt edx,ao<t)= [ [ Py{ao G ds, X„^ G dz)P,{Xt^s G dx) 

Jo J-oo 

and we see that Qy{Xt G •) has a continuous and bounded density function 
given by 

(lt{y,x) = ft{x - y) - ft{y,x), 

where 

ftiy,x)= / Py{aoeds,X„^edz)ft-six- z). 

Jo J-oo 

Then using the Markov property of the meander and the fact that, given 
Et = y, Et+s has law Qy{Xs G •), we see that 

POO 

(19) n*{et£dx,C>t)= n* {et/2 ^ dy,C> t/2)qt,2{y,x) dx 

Jo 

and the existence, continuity and boundedness of pt{-) follows. (This fact is 
also contained in [12], Theorem 6.) Then (17) follows from (13). Of course 
p also has the scaling property (16), so specializing Lemma 6 to the stable 
context and using the identity r(p)r(l — p) = tt/ sinpvr gives (18). □ 

Remark 9. Alternatively, (18) could be deduced from [5], Proposition 
16, page 234. 



We start by showing that we can deduce the behavior of p at zero from 
the fact that it is bounded. 



Proposition 10. (7) holds, viz. p{x) ^Cx°'P as xlO. 

Proof. First we write k{x) for ki(x) with Kt{x) = TT*{t)pt{x) the density 
of n*{et G dx, > t), and note that it suffices to show that 

(20) limx-"''K(x) G (0,oo). 

a; 4-0 

Next, we have seen that Qy{Xt G dx) = Py{Xt G dx,It > 0) = qt{y,x) dx for 
x,y > 0, and by duality (specifically, [5], Theorem 5, page 47), we deduce 
that qt{y,x) =ql{x,y) where 

ql{x, y) dy = P-^{-Xt G dy, 5^ < 0). 
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Using this in the obvious identity 

/•oo 

K2ix) = / K{y)qi{y,x)dy 
Jo 

gives 



K2{x)= I K{y)ql{x,y)dy <c ql{x,y)dy 
Jo Jo 

= cP-^{Si < 0) = cPo{Si < x) - cx'^P. 

(Here, and later, c denotes a generic positive constant whose value can 
change from line to line.) Then, for example, from [8], we recall that the 
law of "—X starting at x > and conditioned to stay positive" has, at time 
1, density 

pI {x,y) ■.= q*i{x,y) — . 

Thus if we write g{y) = K{y)y~°'^ which, by the previous result and scaling, 
is bounded, we see that in the obvious notation. 



K2ix)x-''P = I y~'^PK{y)y'^Pql{x,y)x-"Pdy 

) 

9{y)p*i{x,y)dy 

) 

9{y)Vi {y)dy G (0,cx)), 



where Pi {■) denotes the density of 11 "—X starting at and conditioned 
to stay positive at time 1," and the convergence follows from a result which 
is stated in [9] and proved in [10]. By scaling, we deduce (20) and then (7). 
□ 



Remark 11. Theorem 6 of [12] gives a completely different proof of 
Proposition 10, essentially using random walk approximation. 

Remark 12. Since it is known that the density function of X, start- 
ing from and conditioned to stay positive, is given at time 1 by p^(x) = 
cx°'^^~ k{x) (see, e.g., [8]), we deduce from this, and later from Proposition 
17, that 

(21) p\{x)'^cx°' as 3:4,0 and p\{x) ^ cx~^'^p^^^ asx— >oo. 

In particular we note that p[ and have the same asymptotic behavior at 
0, up to multiplication by a constant. 
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Proposition 13. (6) holds, viz. m{x) ^ Bx'^P"^ as x 4-0. 
Proof. Write = z in (18) and then zx = y to get 

(22) m(x) = ^^^^ r ^^^^^ ax"^-^sinp7r Z'"" p{y) , 

^ ' ^ ' TT {z^-l)P vr y, {y^-x^y 

Since we have seen that p is bounded, 
lim / — < clim 



-5^ A (r - a;")^ - <5W A (r - a;")^ 

-l+<5 

cx^-^^hm / . . =0, 



540 Ji (z" - 1)^ 

as the integral is finite because p < 1. Moreover on ((1 + (5)x,oo), (y° — 
x°')~P < y~°'P{l — (1 + 5)~"^)~P, and, since we know y~°^^p{y) is integrable 
on (0, oo) by Proposition 10, for any (5 > 0, dominated convergence gives 

. . (y^^^^'^y^ ry~'''piy)dy <oc. 
{i+5)x [y — X Jo 

Thus limxiox^~°'fm{x) € (0, oo), and since (4) holds, the limit must be B. 

□ 

Remark 14. It is surprising that the result of Proposition 10 only plays 
a role in our proof of Proposition 13 in the case ap = l. However it is clear 
from (18) that the asymptotic behavior of m and p at infinity are closely 
linked, and we exploit this in the next section. 

4. Proof of the large time results. The first step uses the following ob- 
vious result: 

Lemma 15. Put r = inf{i :et > 1}; then for x> 1 
(23) k{x) 



(24) < 



/ / n*{T edu,er edy)qi_u{y,x) 
Ju=0 Jy=l 

/ / n*{T edu,er €dy)fi_u{x -y). 

Ju=0 Jy=l 



It does not seem possible to deduce the asymptotic behavior of k, from 
(23), but from (24) we can get a useful upper bound. 

Proposition 16. We have 
(25) lim sup2;"'^^K(x) < oo. 
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Proof. Recall that 
(26) lim = A. 

Write the integral in (24) as Ji + and note that 

= /" / {l-u)n*{T edu^Er^dy) 

Ju=0 J\y-x\>x/2,y>l 



A f I {l-u)n*{Tedu,eredy) 

Ju=Q Jy>l 

= A f (1 - u)n*{T e dn) < oo, 

Jn=0 

where we have used (26) and dominated convergence. Also when x is suffi- 
ciently large, 

»1 />3a;/2 



x^+'l2= [ r il-ur^x^+'n*{T€du,eredy)f(j^^) 

Ju=oJx/2 \{^-urJ 

(1 -u)-^x"+^n*(r €du,er£x- dz)f( 



u=0 Jx/2 
1 i-x/2 



u=0 J-x/2 



(1 -li)*? 



But if we put viy) = /J^ J^(^) dz, then for all sufficiently large x, 
x""''"'" sup ?!*(£,-€ X — dzler- = tt^, T = ti) 

0<w<l,-xj2<z<xj2 

a+l V{X- Z-W)dz 

= x sup — <cdz. 

0<w<l-x/2<z<x/2 vil—W) 

So it follows that 

lim sup x"+i/2 <c [ [ (1 - u)-'^n*{T G du)f(—^—) dz 

x^oo Ju=OJ-oo \{^~1''rJ 

1 POO 

/ n*(T G du)f{y) dy = cn*{T < 1) < oo 

u=0 J — oo 

and (25) holds. □ 

To exploit this we specialize (15) to get the following integral equation 
for g: 

X — z 

(27) xg{t,x)= / g{u,z)- ft-u{x - z)dudz. 

Ju=0 Jz=Q t — U 
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We are now in a position to establish the behavior of ■m{x) for large x. 

Proposition 17. If X is any strictly stable process which has positive 
jumps and is not a subordinator, then (5) and (8) hold, viz. m{x) ^ pp{x) ^ 

^^-(a+l) asx— T-OO. 



Proof. Using (17) we can rewrite (27) with t = 1 as 

-1 



(28) 

where, by scaling, 
-1 rx/2 



uP ^pu{z)—^ — ^fi-u{x - z)dudz 



11=0 J z=0 



x{l — u) ' 



J1 + J2 



Jl 



uP-'pu{z){l-- 



1 a+1 



(1-n)^ 



f{{x-z){l-u)-'^)dudz. 



lu=0 J z=0 

Using (26) and dominated convergence gives 

lim Ji = A / uP~^pu{z) dudz = A/ p. 

Ju=0 Jz=0 

Using the scaling property of p and a change of variable, we can write 



J2 = - 



1 


pi px/2 


X 


a+1 , 








pi 


X 


Ju=0 Jy=0 


VTI _ 



x-y 



{l-ur'-'^yf{y{l-u)-^)dudy. 



Now it follows from (25) and the identity, k{x) = 7r*(l)p(x), that the ex- 
pression [:§f]"~^^p{^^^) is bounded above on the range of integration by a 
constant for all sufficiently large x, so that for such x, 



J2<- 



x/2 



uf{l-ur^-'iyf{y{l-ur'i)dudy 
M=0 Jy=Q 

1 rx/2{l-u)'^ 



uP{l-uy-^zf{z)dudz. 

/u=0 J z=0 

Since J^zf{z)dz is 0(l),0(logx), or 0(x^~") according as a > 1, a = 1, or 
a < 1, we check easily that lim^^^oo J2 = 0, and hence that lima;_j.oo x°'^^p{x) - 
A/ p. Putting this into (18) we get 

„,i , , Asmp-n sPds 
lim X m[x) 



p TT 



(1 



A. 



□ 



Remark 18. The reader might like to check that the proof above fails 
if we use the density version of the more obvious identity (14) rather than 
(15). 
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5. Derivatives. It is not difficult to deduce from the known properties of 
/ that p and m are differentiable k times and that these derivatives satisfy 
analogues of (18) and (28); it is also known that 

(29) /('^^ (x) - AfcX-^^+'^+i) as X ^ oo, 

where 

{-)^T{a + l + k)A 



(30) Ak 



r(a + l) 



This suggests that the derivatives of m and p have asymptotic behaviors 
which are consistent with the results of Theorem 1. In fact, we have been 
able to prove such results for large x, but have not been able to settle the 
question for small x. (Details of these results can be supplied by the authors 
on request.) 
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